T-wave inversion through inhomogeneous voltage diffusion within the FK3V cardiac model.

The heart beats are due to the synchronized contraction of cardiomyocytes triggered by a periodic sequence of electrical signals called action potentials, which originate in the sinoatrial node and spread through the heart's electrical system. A large body of work is devoted to modeling the propagation of the action potential and to reproducing reliably its shape and duration. Connection of computational modeling of cells to macroscopic phenomenological curves such as the electrocardiogram has been also intense, due to its clinical importance in analyzing cardiovascular diseases. In this work, we simulate the dynamics of action potential propagation using the three-variable Fenton-Karma model that can account for both normal and damaged cells through a the spatially inhomogeneous voltage diffusion coefficient. We monitor the action potential propagation in the cardiac tissue and calculate the pseudo-electrocardiogram that reproduces the R and T waves. The R-wave amplitude varies according to a double exponential law as a function of the (spatially homogeneous, for an isotropic tissue) diffusion coefficient. The addition of spatial inhomogeneity in the diffusion coefficient by means of a defected region representing damaged cardiac cells may result in T-wave inversion in the calculated pseudo-electrocardiogram. The transition from positive to negative polarity of the T-wave is analyzed as a function of the length and the depth of the defected region.

Enhanced Release Probability without Changes in Synaptic Delay during Analogue-Digital Facilitation.

Neuronal timing with millisecond precision is critical for many brain functions such as sensory perception, learning and memory formation. At the level of the chemical synapse, the synaptic delay is determined by the presynaptic release probability (Pr) and the waveform of the presynaptic action potential (AP). For instance, paired-pulse facilitation or presynaptic long-term potentiation are associated with reductions in the synaptic delay, whereas paired-pulse depression or presynaptic long-term depression are associated with an increased synaptic delay. Parallelly, the AP broadening that results from the inactivation of voltage gated potassium (Kv) channels responsible for the repolarization phase of the AP delays the synaptic response, and the inactivation of sodium (Nav) channels by voltage reduces the synaptic latency. However, whether synaptic delay is modulated during depolarization-induced analogue-digital facilitation (d-ADF), a form of context-dependent synaptic facilitation induced by prolonged depolarization of the presynaptic neuron and mediated by the voltage-inactivation of presynaptic Kv1 channels, remains unclear. We show here that despite Pr being elevated during d-ADF at pyramidal L5-L5 cell synapses, the synaptic delay is surprisingly unchanged. This finding suggests that both Pr- and AP-dependent changes in synaptic delay compensate for each other during d-ADF. We conclude that, in contrast to other short- or long-term modulations of presynaptic release, synaptic timing is not affected during d-ADF because of the opposite interaction of Pr- and AP-dependent modulations of synaptic delay.

A modular organic neuromorphic spiking circuit for retina-inspired sensory coding and neurotransmitter-mediated neural pathways.

Signal communication mechanisms within the human body rely on the transmission and modulation of action potentials. Replicating the interdependent functions of receptors, neurons and synapses with organic artificial neurons and biohybrid synapses is an essential first step towards merging neuromorphic circuits and biological systems, crucial for computing at the biological interface. However, most organic neuromorphic systems are based on simple circuits which exhibit limited adaptability to both external and internal biological cues, and are restricted to emulate only specific the functions of an individual neuron/synapse. Here, we present a modular neuromorphic system which combines organic spiking neurons and biohybrid synapses to replicate a neural pathway. The spiking neuron mimics the sensory coding function of afferent neurons from light stimuli, while the neuromodulatory activity of interneurons is emulated by neurotransmitters-mediated biohybrid synapses. Combining these functions, we create a modular connection between multiple neurons to establish a pre-processing retinal pathway primitive.

A robust balancing mechanism for spiking neural networks.

Dynamical balance of excitation and inhibition is usually invoked to explain the irregular low firing activity observed in the cortex. We propose a robust nonlinear balancing mechanism for a random network of spiking neurons, which works also in the absence of strong external currents. Biologically, the mechanism exploits the plasticity of excitatory-excitatory synapses induced by short-term depression. Mathematically, the nonlinear response of the synaptic activity is the key ingredient responsible for the emergence of a stable balanced regime. Our claim is supported by a simple self-consistent analysis accompanied by extensive simulations performed for increasing network sizes. The observed regime is essentially fluctuation driven and characterized by highly irregular spiking dynamics of all neurons.

Ionic plasmon-polariton in application to neurosignaling.

The diffusive ion current is insufficient to explain the fast saltatory conduction observed in myelinated axons and in pain-sensing C fibers in the human nervous system, where the stimulus signal exhibits a velocity two orders of magnitude greater than the upper limit of ion diffusion velocity, even when the diffusion is accelerated by myelin, as in the discrete cable model including the Hodgkin-Huxley mechanism. The agreement with observations has been achieved in a wave-type model of stimulus signal kinetics via synchronized ion local density oscillations propagating as a wave in axons periodically corrugated by myelin segments in myelinated axons, or by periodically distributed rafts with clusters of Na^{+} channels in C fibers. The resulting so-called plasmon-polariton model for saltatory conduction reveals also the specific role of myelin, which is different from what was previously thought. This can be important for identifying a new target for the future treatment of demyelination diseases.

How neuronal morphology impacts the synchronisation state of neuronal networks.

The biophysical properties of neurons not only affect how information is processed within cells, they can also impact the dynamical states of the network. Specifically, the cellular dynamics of action-potential generation have shown relevance for setting the (de)synchronisation state of the network. The dynamics of tonically spiking neurons typically fall into one of three qualitatively distinct types that arise from distinct mathematical bifurcations of voltage dynamics at the onset of spiking. Accordingly, changes in ion channel composition or even external factors, like temperature, have been demonstrated to switch network behaviour via changes in the spike onset bifurcation and hence its associated dynamical type. A thus far less addressed modulator of neuronal dynamics is cellular morphology. Based on simplified and anatomically realistic mathematical neuron models, we show here that the extent of dendritic arborisation has an influence on the neuronal dynamical spiking type and therefore on the (de)synchronisation state of the network. Specifically, larger dendritic trees prime neuronal dynamics for in-phase-synchronised or splayed-out activity in weakly coupled networks, in contrast to cells with otherwise identical properties yet smaller dendrites. Our biophysical insights hold for generic multicompartmental classes of spiking neuron models (from ball-and-stick-type to anatomically reconstructed models) and establish a connection between neuronal morphology and the susceptibility of neural tissue to synchronisation in health and disease.

Mean-field method for generic conductance-based integrate-and-fire neurons with finite timescales.

The construction of transfer functions in theoretical neuroscience plays an important role in determining the spiking rate behavior of neurons in networks. These functions can be obtained through various fitting methods, but the biological relevance of the parameters is not always clear. However, for stationary inputs, such functions can be obtained without the adjustment of free parameters by using mean-field methods. In this work, we expand current Fokker-Planck approaches to account for the concurrent influence of colored and multiplicative noise terms on generic conductance-based integrate-and-fire neurons. We reduce the resulting stochastic system through the application of the diffusion approximation to a one-dimensional Langevin equation. An effective Fokker-Planck is then constructed using Fox Theory, which is solved numerically using a newly developed double integration procedure to obtain the transfer function and the membrane potential distribution. The solution is capable of reproducing the transfer function and the stationary voltage distribution of simulated neurons across a wide range of parameters. The method can also be easily extended to account for different sources of noise with various multiplicative terms, and it can be used in other types of problems in principle.

Intracellular ion accumulation in the genesis of complex action potential dynamics under cardiac diseases.

Intracellular ions, including sodium (Na^{+}), calcium (Ca^{2+}), and potassium (K^{+}), etc., accumulate slowly after a change of the state of the heart, such as a change of the heart rate. The goal of this study is to understand the roles of slow ion accumulation in the genesis of cardiac memory and complex action-potential duration (APD) dynamics that can lead to lethal cardiac arrhythmias. We carry out numerical simulations of a detailed action potential model of ventricular myocytes under normal and diseased conditions, which exhibit memory effects and complex APD dynamics. We develop a low-dimensional iterated map (IM) model to describe the dynamics of Na^{+}, Ca^{2+}, and APD and use it to uncover the underlying dynamical mechanisms. The development of the IM model is informed by simulation results under the normal condition. We then use the IM model to perform linear stability analyses and computer simulations to investigate the bifurcations and complex APD dynamics, which depend on the feedback loops between APD and intracellular Ca^{2+} and Na^{+} concentrations and the steepness of the APD response to the ion concentrations. When the feedback between APD and Ca^{2+} concentration is positive, a Hopf bifurcation leading to periodic oscillatory behavior occurs as the steepness of the APD response to the ion concentrations increases. The negative feedback loop between APD and Na^{+} concentration is required for the Hopf bifurcation. When the feedback between APD and Ca^{2+} concentration is negative, period-doubling bifurcations leading to high periodicity and chaos occurs. In this case, Na^{+} accumulation plays little role in the dynamics. Finally, we carry out simulations of the detailed action potential model under two diseased conditions, which exhibit steep APD responses to ion concentrations. Under both conditions, Hopf bifurcations leading to slow oscillations or period-doubling bifurcations leading to high periodicity and chaotic APD dynamics occur, depending on the strength of the ion pump-Na^{+}-Ca^{2+} exchanger. Using functions reconstructed from the simulation data, the IM model accurately captures the bifurcations and dynamics under the two diseased conditions. In conclusion, besides using computer simulations of a detailed high-dimensional action-potential model to investigate the effects of slow ion accumulation and short-term memory on bifurcations and genesis of complex APD dynamics in cardiac myocytes under diseased conditions, this study also provides a low-dimensional mathematical tool, i.e., the IM model, to allow stability analyses for uncovering the underlying mechanisms.

Linear and nonlinear integrate-and-fire neurons driven by synaptic shot noise with reversal potentials.

The steady-state firing rate and firing-rate response of the leaky and exponential integrate-and-fire models receiving synaptic shot noise with excitatory and inhibitory reversal potentials is examined. For the particular case where the underlying synaptic conductances are exponentially distributed, it is shown that the master equation for a population of such model neurons can be reduced from an integrodifferential form to a more tractable set of three differential equations. The system is nevertheless more challenging analytically than for current-based synapses: where possible, analytical results are provided with an efficient numerical scheme and code provided for other quantities. The increased tractability of the framework developed supports an ongoing critical comparison between models in which synapses are treated with and without reversal potentials, such as recently in the context of networks with balanced excitatory and inhibitory conductances.

Computational models of compound nerve action potentials: Efficient filter-based methods to quantify effects of tissue conductivities, conduction distance, and nerve fiber parameters.

BACKGROUND: Peripheral nerve recordings can enhance the efficacy of neurostimulation therapies by providing a feedback signal to adjust stimulation settings for greater efficacy or reduced side effects. Computational models can accelerate the development of interfaces with high signal-to-noise ratio and selective recording. However, validation and tuning of model outputs against in vivo recordings remains computationally prohibitive due to the large number of fibers in a nerve. METHODS: We designed and implemented highly efficient modeling methods for simulating electrically evoked compound nerve action potential (CNAP) signals. The method simulated a subset of fiber diameters present in the nerve using NEURON, interpolated action potential templates across fiber diameters, and filtered the templates with a weighting function derived from fiber-specific conduction velocity and electromagnetic reciprocity outputs of a volume conductor model. We applied the methods to simulate CNAPs from rat cervical vagus nerve. RESULTS: Brute force simulation of a rat vagal CNAP with all 1,759 myelinated and 13,283 unmyelinated fibers in NEURON required 286 and 15,860 CPU hours, respectively, while filtering interpolated templates required 30 and 38 seconds on a desktop computer while maintaining accuracy. Modeled CNAP amplitude could vary by over two orders of magnitude depending on tissue conductivities and cuff opening within experimentally relevant ranges. Conduction distance and fiber diameter distribution also strongly influenced the modeled CNAP amplitude, shape, and latency. Modeled and in vivo signals had comparable shape, amplitude, and latency for myelinated fibers but not for unmyelinated fibers. CONCLUSIONS: Highly efficient methods of modeling neural recordings quantified the large impact that tissue properties, conduction distance, and nerve fiber parameters have on CNAPs. These methods expand the computational accessibility of neural recording models, enable efficient model tuning for validation, and facilitate the design of novel recording interfaces for neurostimulation feedback and understanding physiological systems.

Synchrony in Networks of Type 2 Interneurons Is More Robust to Noise with Hyperpolarizing Inhibition Compared to Shunting Inhibition in Both the Stochastic Population Oscillator and the Coupled Oscillator Regimes.

Synchronization in the gamma band (25-150 Hz) is mediated by PV+ inhibitory interneurons, and evidence is accumulating for the essential role of gamma oscillations in cognition. Oscillations can arise in inhibitory networks via synaptic interactions between individual oscillatory neurons (mean-driven) or via strong recurrent inhibition that destabilizes the stationary background firing rate in the fluctuation-driven balanced state, causing an oscillation in the population firing rate. Previous theoretical work focused on model neurons with Hodgkin's Type 1 excitability (integrators) connected by current-based synapses. Here we show that networks comprised of simple Type 2 oscillators (resonators) exhibit a supercritical Hopf bifurcation between synchrony and asynchrony and a gradual transition via cycle skipping from coupled oscillators to stochastic population oscillator (SPO), as previously shown for Type 1. We extended our analysis to homogeneous networks with conductance rather than current based synapses and found that networks with hyperpolarizing inhibitory synapses were more robust to noise than those with shunting synapses, both in the coupled oscillator and SPO regime. Assuming that reversal potentials are uniformly distributed between shunting and hyperpolarized values, as observed in one experimental study, converting synapses to purely hyperpolarizing favored synchrony in all cases, whereas conversion to purely shunting synapses made synchrony less robust except at very high conductance strengths. In mature neurons the synaptic reversal potential is controlled by chloride cotransporters that control the intracellular concentrations of chloride and bicarbonate ions, suggesting these transporters as a potential therapeutic target to enhance gamma synchrony and cognition.

An operating principle of the cerebral cortex, and a cellular mechanism for attentional trial-and-error pattern learning and useful classification extraction.

A feature of the brains of intelligent animals is the ability to learn to respond to an ensemble of active neuronal inputs with a behaviorally appropriate ensemble of active neuronal outputs. Previously, a hypothesis was proposed on how this mechanism is implemented at the cellular level within the neocortical pyramidal neuron: the apical tuft or perisomatic inputs initiate "guess" neuron firings, while the basal dendrites identify input patterns based on excited synaptic clusters, with the cluster excitation strength adjusted based on reward feedback. This simple mechanism allows neurons to learn to classify their inputs in a surprisingly intelligent manner. Here, we revise and extend this hypothesis. We modify synaptic plasticity rules to align with behavioral time scale synaptic plasticity (BTSP) observed in hippocampal area CA1, making the framework more biophysically and behaviorally plausible. The neurons for the guess firings are selected in a voluntary manner via feedback connections to apical tufts in the neocortical layer 1, leading to dendritic Ca2+ spikes with burst firing, which are postulated to be neural correlates of attentional, aware processing. Once learned, the neuronal input classification is executed without voluntary or conscious control, enabling hierarchical incremental learning of classifications that is effective in our inherently classifiable world. In addition to voluntary, we propose that pyramidal neuron burst firing can be involuntary, also initiated via apical tuft inputs, drawing attention toward important cues such as novelty and noxious stimuli. We classify the excitations of neocortical pyramidal neurons into four categories based on their excitation pathway: attentional versus automatic and voluntary/acquired versus involuntary. Additionally, we hypothesize that dendrites within pyramidal neuron minicolumn bundles are coupled via depolarization cross-induction, enabling minicolumn functions such as the creation of powerful hierarchical "hyperneurons" and the internal representation of the external world. We suggest building blocks to extend the microcircuit theory to network-level processing, which, interestingly, yields variants resembling the artificial neural networks currently in use. On a more speculative note, we conjecture that principles of intelligence in universes governed by certain types of physical laws might resemble ours.

Silencing CA1 pyramidal cells output reveals the role of feedback inhibition in hippocampal oscillations.

The precise temporal coordination of neural activity is crucial for brain function. In the hippocampus, this precision is reflected in the oscillatory rhythms observed in CA1. While it is known that a balance between excitatory and inhibitory activity is necessary to generate and maintain these oscillations, the differential contribution of feedforward and feedback inhibition remains ambiguous. Here we use conditional genetics to chronically silence CA1 pyramidal cell transmission, ablating the ability of these neurons to recruit feedback inhibition in the local circuit, while recording physiological activity in mice. We find that this intervention leads to local pathophysiological events, with ripple amplitude and intrinsic frequency becoming significantly larger and spatially triggered local population spikes locked to the trough of the theta oscillation appearing during movement. These phenotypes demonstrate that feedback inhibition is crucial in maintaining local sparsity of activation and reveal the key role of lateral inhibition in CA1 in shaping circuit function.

Local origin of excitatory-inhibitory tuning equivalence in a cortical network.

The interplay between excitation and inhibition determines the fidelity of cortical representations. The receptive fields of excitatory neurons are often finely tuned to encoded features, but the principles governing the tuning of inhibitory neurons remain elusive. In this study, we recorded populations of neurons in the mouse postsubiculum (PoSub), where the majority of excitatory neurons are head-direction (HD) cells. We show that the tuning of fast-spiking (FS) cells, the largest class of cortical inhibitory neurons, was broad and frequently radially symmetrical. By decomposing tuning curves using the Fourier transform, we identified an equivalence in tuning between PoSub-FS and PoSub-HD cell populations. Furthermore, recordings, optogenetic manipulations of upstream thalamic populations and computational modeling provide evidence that the tuning of PoSub-FS cells has a local origin. These findings support the notion that the equivalence of neuronal tuning between excitatory and inhibitory cell populations is an intrinsic property of local cortical networks.

Cortical cell assemblies and their underlying connectivity: An in silico study.

Recent developments in experimental techniques have enabled simultaneous recordings from thousands of neurons, enabling the study of functional cell assemblies. However, determining the patterns of synaptic connectivity giving rise to these assemblies remains challenging. To address this, we developed a complementary, simulation-based approach, using a detailed, large-scale cortical network model. Using a combination of established methods we detected functional cell assemblies from the stimulus-evoked spiking activity of 186,665 neurons. We studied how the structure of synaptic connectivity underlies assembly composition, quantifying the effects of thalamic innervation, recurrent connectivity, and the spatial arrangement of synapses on dendrites. We determined that these features reduce up to 30%, 22%, and 10% of the uncertainty of a neuron belonging to an assembly. The detected assemblies were activated in a stimulus-specific sequence and were grouped based on their position in the sequence. We found that the different groups were affected to different degrees by the structural features we considered. Additionally, connectivity was more predictive of assembly membership if its direction aligned with the temporal order of assembly activation, if it originated from strongly interconnected populations, and if synapses clustered on dendritic branches. In summary, reversing Hebb's postulate, we showed how cells that are wired together, fire together, quantifying how connectivity patterns interact to shape the emergence of assemblies. This includes a qualitative aspect of connectivity: not just the amount, but also the local structure matters; from the subcellular level in the form of dendritic clustering to the presence of specific network motifs.

Beyond a Transmission Cable-New Technologies to Reveal the Richness in Axonal Electrophysiology.

The axon is a neuronal structure capable of processing, encoding, and transmitting information. This assessment contrasts with a limiting, but deeply rooted, perspective where the axon functions solely as a transmission cable of somatodendritic activity, sending signals in the form of stereotypical action potentials. This perspective arose, at least partially, because of the technical difficulties in probing axons: their extreme length-to-diameter ratio and intricate growth paths preclude the study of their dynamics through traditional techniques. Recent findings are challenging this view and revealing a much larger repertoire of axonal computations. Axons display complex signaling processes and structure-function relationships, which can be modulated via diverse activity-dependent mechanisms. Additionally, axons can exhibit patterns of activity that are dramatically different from those of their corresponding soma. Not surprisingly, many of these recent discoveries have been driven by novel technology developments, which allow for in vitro axon electrophysiology with unprecedented spatiotemporal resolution and signal-to-noise ratio. In this review, we outline the state-of-the-art in vitro toolset for axonal electrophysiology and summarize the recent discoveries in axon function it has enabled. We also review the increasing repertoire of microtechnologies for controlling axon guidance which, in combination with the available cutting-edge electrophysiology and imaging approaches, have the potential for more controlled and high-throughput in vitro studies. We anticipate that a larger adoption of these new technologies by the neuroscience community will drive a new era of experimental opportunities in the study of axon physiology and consequently, neuronal function.

Impact on backpropagation of the spatial heterogeneity of sodium channel kinetics in the axon initial segment.

In a variety of neurons, action potentials (APs) initiate at the proximal axon, within a region called the axon initial segment (AIS), which has a high density of voltage-gated sodium channels (NaVs) on its membrane. In pyramidal neurons, the proximal AIS has been reported to exhibit a higher proportion of NaVs with gating properties that are "right-shifted" to more depolarized voltages, compared to the distal AIS. Further, recent experiments have revealed that as neurons develop, the spatial distribution of NaV subtypes along the AIS can change substantially, suggesting that neurons tune their excitability by modifying said distribution. When neurons are stimulated axonally, computational modelling has shown that this spatial separation of gating properties in the AIS enhances the backpropagation of APs into the dendrites. In contrast, in the more natural scenario of somatic stimulation, our simulations show that the same distribution can impede backpropagation, suggesting that the choice of orthodromic versus antidromic stimulation can bias or even invert experimental findings regarding the role of NaV subtypes in the AIS. We implemented a range of hypothetical NaV distributions in the AIS of three multicompartmental pyramidal cell models and investigated the precise kinetic mechanisms underlying such effects, as the spatial distribution of NaV subtypes is varied. With axonal stimulation, proximal NaV availability dominates, such that concentrating right-shifted NaVs in the proximal AIS promotes backpropagation. However, with somatic stimulation, the models are insensitive to availability kinetics. Instead, the higher activation threshold of right-shifted NaVs in the AIS impedes backpropagation. Therefore, recently observed developmental changes to the spatial separation and relative proportions of NaV1.2 and NaV1.6 in the AIS differentially impact activation and availability. The observed effects on backpropagation, and potentially learning via its putative role in synaptic plasticity (e.g. through spike-timing-dependent plasticity), are opposite for orthodromic versus antidromic stimulation, which should inform hypotheses about the impact of the developmentally regulated subcellular localization of these NaV subtypes.

Fractional order memcapacitive neuromorphic elements reproduce and predict neuronal function.

There is an increasing need to implement neuromorphic systems that are both energetically and computationally efficient. There is also great interest in using electric elements with memory, memelements, that can implement complex neuronal functions intrinsically. A feature not widely incorporated in neuromorphic systems is history-dependent action potential time adaptation which is widely seen in real cells. Previous theoretical work shows that power-law history dependent spike time adaptation, seen in several brain areas and species, can be modeled with fractional order differential equations. Here, we show that fractional order spiking neurons can be implemented using super-capacitors. The super-capacitors have fractional order derivative and memcapacitive properties. We implemented two circuits, a leaky integrate and fire and a Hodgkin-Huxley. Both circuits show power-law spiking time adaptation and optimal coding properties. The spiking dynamics reproduced previously published computer simulations. However, the fractional order Hodgkin-Huxley circuit showed novel dynamics consistent with criticality. We compared the responses of this circuit to recordings from neurons in the weakly-electric fish that have previously been shown to perform fractional order differentiation of their sensory input. The criticality seen in the circuit was confirmed in spontaneous recordings in the live fish. Furthermore, the circuit also predicted long-lasting stimulation that was also corroborated experimentally. Our work shows that fractional order memcapacitors provide intrinsic memory dependence that could allow implementation of computationally efficient neuromorphic devices. Memcapacitors are static elements that consume less energy than the most widely studied memristors, thus allowing the realization of energetically efficient neuromorphic devices.

Adaptive myelination causes slow oscillations in recurrent neural loops.

The brain is known to be plastic, i.e., capable of changing and reorganizing as it develops and accumulates experience. Recently, a novel form of brain plasticity was described which is activity-dependent myelination of nerve fibers. Since the speed of propagation of action potentials along axons depends significantly on their degree of myelination, this process leads to adaptive change of axonal delays depending on the neural activity. To understand the possible influence of the adaptive delays on the behavior of neural networks, we consider a simple setup, a neuronal oscillator with delayed feedback. We show that introducing the delay plasticity into this circuit can lead to the occurrence of slow oscillations which are impossible with a constant delay.